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and midrapidity pt distributions of hadrons and direct photons 
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In the frame of 2+1 Bjorken hydrodynamics we describe 
simuhaneously hadron yields and midrapidity pt distributions 
of direct photons, unflavored, and strange hadrons measured 
in Pb+Pb collisions at SpS. We find, that a reasonably quanti- 
tative description of these data is possible only if we introduce 
some radial velocity at the beginning of the one-fluid hydro- 
dynamic stage. We fix uniquely all parameters of this model, 
and estimate an initial temperature of the just thermalized 
state Tin = 200^20 -^sV in the case of EoS with phase tran- 
sition to Quark-Gluon Plasma, and Tin ~ 230 [j^?? MeV in the 
pure hadronic scenario. We conclude that SpS data consid- 
ered in our study can not distinguish between a production 
of QGP, and its absence. 

A finite volume of rapidly expanding Quark-Gluon 
Plasma (QGP) can be possibly created for a very short 
time in ultrarelativistic heavy ion collisions as an almost 
locally equihbrated thermodynamic state of deconfined 
quarks and gluons with strong residual interaction. Tak- 
ing into account a very short time available for QGP 
formation in this process a necessary condition of its cre- 
ations is that an initial (thermalization) temperature Tin 
have to be visibly higher than the quark-hadron phase 
transition temperature, which is estimated from lattice 
calculations to be T^ = 150- 170A/ey [1]. Therefore an 
estimation of a temperature at which the local thermal 
equilibrium is just established - thermalization temper- 
ature - is of great importance. 

There are a lot of attempts to extract this tempera- 
ture from direct photon momentum distributions mea- 
sured in S+Au [2] and Pb+Pb [3] collisions. However, 
to extract the thermalization temperature from experi- 
mental data reasonably well, it is necessary to separate 
accurately the collective and chaotic motions. To do this 
one needs to know the evolution of the collision start- 
ing from the first interactions till the final particle stage. 
Unfortunately, the most part of the estimations of the 
thermalization temperature [4,5] are based on hydrody- 
namic description with too simplified initial conditions 
(zero initial radial collective velocity), and other model 
parameters fixed partially by some experimental data, 
and partially - ad hoc from more or less reasonable con- 
siderations. Only in the recent paper [6] a first attempt 
to estimate the influence of nonzero initial radial velocity 
was undertaken. 

In contrast to this, in the present letter we use initial 
conditions for hydrodynamic model in the most general 
form, describe simultaneously a large amount of avail- 
able experimental data, measured in Pb+Pb collisions at 



SpS (hadron yields, midrapidity pt distributions of direct 
photons, unflavored and strange hadrons), fix uniquely 
all model parameters from these data and obtain reli- 
able estimate of the thermalization temperature with its 
uncertainty concerned with experimental data bars. 

In this way we find, that if we use a realistic equa- 
tion of state (EoS) of hadronic gas incorporating con- 
tributions of all known hadrons and resonances, then to 
describe midrapidity pt distributions of final hadrons it 
is necessary to assume an essentially nonzero initial ra- 
dial collective velocity at the beginning of the one-fluid 
hydrodynamic stage. This velocity is created during the 
pre-equilibrium (two-fluid) stage of the collision due to 
significant radial energy density gradient coming from the 
shapes of the colliding nuclei. Nevertheless, in all previ- 
ous estimations of the thermalization temperature this 
initial radial velocity was absolutely neglected, and only 
the simplest case (zero initial velocity) was considered. 

Evaluating the direct photon yield we assume the fol- 
lowing scenario of the collision. Initially, on the pre-equi- 
librium stage, colliding nuclei penetrate through each 
other like two fluids with a friction. On this stage the 
'prompt' photons are emitted and radial collective veloc- 
ity is created. If the thermalization time is not negligibly 
small then the both these effects should be taken into 
account. The next stage is a one- fluid- like expansion of 
the locally thermalized matter with an emission of the 
'thermal' photons and its decay onto final hadrons. 

To begin with we calculate contribution of the prompt 
photons. We take into account Compton, annihilation 
and bremsstrahlung processes. We use GRV-94 structure 
functions [7], two-loop expression for as and K-factor, 
K=2, accounting higher order corrections. Comparing 
calculated yield of prompt photons in pp and pA colli- 
sions with experimental data at ^/s = 19 GeV - the clos- 
est available cm. energy to Pb+Pb at SpS (-^/s — 17.3 
GeV) we find, that both experimental data and theoret- 
ical predictions have approximately the same slope, but 
difference in normalization reaches factor « 7. Nowadays 
this difference is attributed to the intrinsic transverse mo- 
mentum of colliding partons [8] . In this letter we account 
these corrections by introducing a proper factor, assum- 
ing that it is independent on ^/s while moving from pA to 
Pb+Pb energy. To go from pA to AA collisions we need 
the number of nucleon-nuclcon collisions occurred during 
the interpenetration of the colliding nuclei. In our calcu- 
lations we use Glauber approximation, what gives in the 
case of Pb+Pb collision Nqi ^ 3 ■ A. Comparison with 



WA98 data [3] shows, that prompt photons contribute 
approximately 30% of total direct photon yield. 

To evaluate yield of thermal photons we integrate aver- 
age emission rate of direct photons from unit volume per 
unit time over a space-time evolution of hot matter calcu- 
lated in a one-fluid hydrodynamics. In the present anal- 
ysis we use recently calculated emission rates account- 
ing next to leading effects in QGP [9] and dominating 
contribution in hadronic gas - reaction np ^ ai ^ t:"/ 
[10]. For description of the evolution we use 2+1 Bjorkcn 
hydrodynamics, which works sufficiently well at midra- 
pidity It has several parameters which should be fixed 
to define the evolution: a) initial conditions: initial tem- 
perature Tin, initial time Tin, initial energy density and 
velocity distributions; b) EoS of hot matter: transition 
temperature Tc (if EoS includes phase transition) , chem- 
ical freeze-out temperature Teh and thermal freeze-out 
temperature Tf. 

For the shape of energy density distribution we use 
Woods-Saxon distribution, integrated along beam axis. 
As for initial distribution of radial velocity, it deserve 
separate discussion. As we find, to describe momentum 
distributions of final hadrons using realistic EoS it is nec- 
essary to introduce some radial velocity already at the 
beginning of the onc-fiuid hydrodynamics. This velocity 
is created during interpenetration of colliding nuclei due 
to significant radial energy density gradient, originated 
from spherical shapes of colliding nuclei. Two-fluid hy- 
drodynamics or kinetic models might predict value and 
radial dependence of this collective velocity but, in this 
letter we use 'phemcnologic' approach and use a simple 
linear distribution v{r) — (r/R) ■Q{R—r)- V"^°'^ charac- 
terized by one new parameter - V™""^ - the velocity on 
the surface of the just thermalized system. We find, that 
our results are not sensitive to reasonable variations of 
the shape of this distribution. 

Concerning equation of state we consider two differ- 
ent cases: EoS with phase transition to QGP, and pure 
hadronic EoS. In the case of phase transition for QGP 
phase we use EoS of ideal gas of massless quarks and 
gluons with degeneracy g — 41.5, what corresponds to 
QGP consisting of 2.5 massless quark flavors and gluons. 
To calculate EoS of hadronic phase both in the case of 
phase transition and pure hadronic scenario we take into 
account contributions of all hadrons listed in the Parti- 
cle Data Book [11]. Inclusion of heavy hadrons leads to 
significant reduction of speed of sound in the hadronic 
phase: v"^ ~ 0.15 in our case should be compared with 
Wg « 0.33 for massless hadronic gas or v"^ « 0.25 for 
pure pion gas. The decreasing of speed of sound results 
in decreasing of acceleration in hadronic gas and stepper 
slopes oi pt distributions of final hadrons [12]. 

To model freeze-out of hadronic gas we use two differ- 
ent parameters: chemical (Teh) and thermal {Tf) freeze- 
out temperatures. It is known, that due to large dif- 
ference in cross-sections of elastic and inelastic reactions 



in hadronic gas, the chemical freeze-out takes place be- 
fore the thermal one. Following paper [13] we use the 
following values of the temperature of chemical freeze- 
out, and baryon, strange and /3 chemical potentials - 
Teh == 168 MeV, ^ib = 266 MeV, p.s = TlMeV and 
/i3 =: —5 MeV. As a result we well describe yields of 
TT, K, ij, p. A, S and H, hadrons measured by NA44 [14], 
NA49 [15], WA97 [16] and WA98 [17] coUaborations in 
Pb+Pb colfisions at SpS energy. 

A description of our results we start from pure hadronic 
scenario. Despite the belief that hadronic gas can not ex- 
ist at temperatures higher than ^ 200 MeV (because of 
the essential overlapping of hadrons in it) we adopt its 
existence and consider the consequences of this assump- 
tion from the point of view of description of direct photon 
production. To begin with, we assume zero initial radial 
velocity and try to describe momentum distributions of 
final hadrons. We choose some initial temperature and 
calculate initial time to describe multiplicity of final 7r°. 
After that wc perform the hydrodynamic calculations 
with various values of chemical and thermal freeze-out 
temperatures and obtain xj of the fit of the calculated 
to the experimental tt*^ distributions as a function of Teh 
and Tf. We find, that it is impossible to describe midra- 
pidity pt distributions of tt" within any reasonable set 
of model parameters: because of the much steeper slope 
of the calculated pt distribution of tt" with respect to 
the experimental one, we obtain xi ^ 10^ /point. We 
find the same situation for all initial temperatures in the 
range 160 MeV < Tin < 300 MeV. A reasonable way to 
improve this situation is to introduce some initial radial 
velocity to the one-fluid hydrodynamic stage. 
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FIG. 1. Levels of constant Xvr of the fit of the ratio of 
calculated to measured pt distributions of tt" as a function of 
maximal initial velocity (j/™"^) and thermal freeze-out tem- 
peratures (T/) for pure hadronic gas EoS, Tin = 230 MeV. 



On the next step, we repeat calculations with nonzero 
initial radial velocity. Again we choose initial tempera- 
ture, fix initial time from the multiplicity considerations, 
fix temperature of chemical freeze-out Teh = 168 MeV 
to describe hadron yields and plot levels of constant xj 
in the Tf — \/™°^ coordinates - see fig. 1. The function 
xi{Tf,V"^'^^) has a shape of a valley with steep walls, 
while any point on its bottom (thick line), corresponds 
to a good description of midrapidity pt distributions of 
TT^ with xi ~ i/point. To fix a point on this line we use 
spectra of heavier hadrons {K, p, A, S). It was shown 
[18] that spectra of the heavier hadrons can be described 
well if we assume freeze-out temperature Tf = 120 MeV. 

The description of the hadron pt distributions at 
midrapidity obtained in this way is shown on the fig. 2. 
One can see, that shapes of all available experimental 
spectra: tt" [17], K and A [19], protons [20], and S [21] 
are described very well, although we overestimate total 
yields of p and A in this rough model. 




FIG. 2. Comparison of experimental pt distributions 
of several hadrons with hydrodynamic predictions for pure 
hadronic gas EoS with Tin = 230 MeV and initial radial ve- 
locity V'""'' = 0.26. 

To fix finally all model parameters we estimate Tm 
using Pt distribution of direct photons. We depict on 
the fig. 3 the sum of the pt distribution of prompt and 
thermal photons, calculated at three initial temperatures. 
For each initial temperature we repeat calculations, de- 
scribed above and obtain the following sets of parameters: 
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We find, an agreement with experimental data at ini- 
tial temperature T,„ = 230+!^° Me^: T,„ = 230 MeV 
corresponds to xi — 3.7/point, Xy — 0.&7 /point, 
while curves, calculated for Tin — 200 MeV (x^ — 
4.9/point, xi = 0.67 /point), and r,„ = 300 MeV (xl = 
3.5/point, Xj — 2.3/point) are still within experimental 
errors. The reason of this weak sensitivity to the initial 
temperature is following. As far as we introduce the ini- 
tial radial velocity, all time stages contribute with more 
or less the same effective slopes. But, the difference be- 
tween initial time for Tin = 300 MeV - Tin = 0.12 fm/c, 
and for Tin = 230 MeV - Tin = 0.6 fm/c is not large, and 
therefore contribution of this stage is not very important 
with respect to the subsequent evolution. 




FIG. 3. Yield of direct photons in Pb-|-Pb collisions at 
SpS energy with pure hadronic gas EoS for different Ti„: 200 
(solid line), 230 (dashed line) and 300 MeV (dotted line). 
Dots - WA98 data. 

We confirm the result found in [5] - 'reach' hadron 
gas EoS allows to describe the experimental pt distribu- 
tion of direct photons [3] within pure hadronic scenario. 
Our best fit value of Tin — 230 ^gg MeV is smaller, than 
the initial temperature estimated in [5] Tin ~ 260 MeV. 
There are two reasons of this discrepancy: first is, the 
initial radial velocity which we included to describe final 
hadron spectra. In contrast to our approach, authors of 
[5] made not attempt to describe hadron spectra, and did 
not consider a nonzero initial radial velocity. The second 
reason is different normalization to final hadron multi- 
plicity. The authors of [5] used the well known Bjorken 
formula, relating initial time and temperature with final 
multiplicity (formula (5) in [5]). But, this formula is valid 
only for massless final particles and being applied to real 
hadronic gas leads to an underestimation of the initial 



time and, as a consequence, of the thermal photon yield. 

Now let us consider EoS incorporating first order phase 
transition. In this case there is an additional parameter 
Tc which value is estimated from lattice QCD simulations 
to be 150 - 170 MeV [1]. On the other hand Tc > T^h- 
Therefore we use Tc ~ Teh = 168 MeV to consider an 
extreme case with the most intensive QGP production. 

As in the case of pure hadronic gas, first we try to de- 
scribe experimental pt distribution of tt" with zero initial 
radial velocity, and find almost the same situation: for 
any reasonable set of model parameters the calculated pt 
distribution of n^ goes much steeper, and result in unrea- 
sonably large xi ~ 10^ /point. So, we find again, that it 
is necessary to introduce initial radial velocity to describe 
momentum distributions of final hadrons. Repeating the 
same procedure as for pure hadronic gas EoS we obtain 
the following sets of model parameters: 
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We find the best agreement with experimental data at 



T, 



200 



+40 



20 MeV (xj = 3.3/pomi, x' = 0.8/point), 
what means very short QGP phase and much more pro- 
longed mixed phase. 




FIG. 4. Direct photon yields calculated for EoS with phase 
transition and linear initial velocity distribution at Tin ~ 180 
MeV (solid line), T„ = 200 MeV (dashed line), T„ = 230 
MeV (dotted line) and Tn = 250 MeV (dash-dotted line). 

To conclude, we use 2+1 Bjorken hydrodynamics to 
analyze hadron yields, and midrapidity pt distributions 



of hadrons and direct photons, measured in Pb+Pb col- 
lision at SpS energy. We find, that calculations with zero 
initial velocity, and EoS of hadronic gas accounting con- 
tribution of all known hadrons result in too soft pt distri- 
butions of final hadrons, and a reasonable way to describe 
experimental slopes is to introduce nonzero radial veloc- 
ity at the beginning of the one-fiuid hydrodynamic stage. 
This initial radial velocity, which might be produced dur- 
ing the pre-equilibrium stage because of the strong radial 
energy density gradient originated from spherical shapes 
of colliding nuclei, is not small (« 0.3 c near the sur- 
face) . Therefore the pre-equilibrium stage is expected to 
be long enough, and, in particular the 'prompt photons' 
have to be taken into account. 

An introduction of the initial radial velocity to this 
model allows to describe simultaneously direct photon 
and hadron pt distributions for the both EoS includ- 
ing phase transition and without it. We estimate an 
initial temperature of the just thcrmalized state Tin = 
200^20 MeV in the case of EoS with phase transition to 
QGP, and T,„ = 230 +3Q MeV in the pure hadronic sce- 
nario. Because the lower bounds of the estimated ther- 
malization temperature are very close to the expected 
transition temperature we conclude that the SpS data 
considered in this letter can not distinguish between a 
production of QGP, and its absence: The initial temper- 
ature in central Pb-|-Pb collisions at SpS is too small, and 
experimental errors for the direct photon pt distributions 
are too large. 
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